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^ INTRODUCTION 

m 

O Functional integral methods are one way of discussing the physics of interacting fermions 

^ which in many cases turns out to be particularly transparent and appealing. In particular in 

O cases with a broken symmetry the use of the Hubbard-Stratonovich transformation^ often 

^ allows one to formulate the problem in a way that is both physically transparent and systematic. 

I here discuss some applications to the Hubbard model, both attractive and repulsive. 

^ The standard Hubbard Hamiltonian has the form 

s 

^ H = -tY, {O'ta'^r'a + ^-C.) + t/ ^ rhr^rirl , (1) 

^ {rr') r 

o 

O where a],^ creates a fermion at site r with spin projection a,tis the nearest neighbor hopping 

integral, U is the onsite interaction (either attractive or repulsive) , and (rr') indicates 
summation over nearest-neighbor bonds, each bond being counted once. Introducing a 
spinor notation via 



the partition function can be written as a functional integral over Grassmann variables: 

Z = j VAe-^"-^'"' . (3) 
The free and interaction parts of the action are respectively 

So= [^dr\Y^ %(dr - ^i)^r - t E + c.c.) \ , (4) 

^° I r irr>) J 



and 



r0 

Sint = U drY.i^t^'^ri'^rl'^r^ ■ (5) 



In the following, the interaction term will be treated using a Hubbard-Stratonovich decom- 
position, which however is defined slightly differently for the attractive and repulsive cases. 



ATTRACTIVE HUBBARD MODEL 

In this case we use the identity 



(6) 



where A is a complex variable, and e is the "thickness" of the time slices used to define 
the functional integral.^ Inserting this at each point in space and at each of the time slices 
the quartic term Sint now has been decomposed into a purely bilinear form, but of couse 
there is now a functional integration over A(r, r). In the saddle point approximation for A 
one recovers BCS theory, and in the following I will derive the effective action for the low- 
energy excitations around the BCS ground state. For the moment I will consider a "neutral 
superconductor" (i.e. a superfluid) and therefore neglect the coupling to the electromagnetic 
filed. The standard way to do the expansion is to write 

A(r,T) = Ao + 5(r,T) . (7) 

where Aq is the BCS gap parameter, and then to expand in 6. However, in this way amplitude 
and phase excitations are mixed up and some effort is needed to disentangle both and to obtain 
explicitly the important phase excitations. The derivation can be considerably simplified by 
performing first a gauge transformation on the fermions: one writes 

A(r,T) = |A(r,r)|e^^(^'^^ , (8) 

and defines new fermion fields $ by 

= e*^(^'^)/2$r . (9) 

The action then takes the form 

S = Sbcs{^, Ao) + ^) + Se{^, | A| - Aq) . (10) 

Here Sbcs is the mean-field action, contains the amplitude fluctuations and is of minor 
importance as long as one is well below Tc and in the weak-coupling limit (so that the 
coherence length is large). The important term is the contribution from the phase degrees of 
freedom: 



■^0 \ ^ r (r,r') 



(11) 



This is now explicitly of at least first order in the derivatives of (p, and consequently in order 
to obtain the lowest order effective action, one has to expand the action to second order in 



S^p. The result is 



S^fii^) = \ [ dr I d'riK^' + psi^if)'} . (12) 
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The coefficient k is related to the compressibility of the fermions in the BCS state and given 
by 

K = \ {nn)q^^^Q , (13) 
and the superfluid density is 

Ps = Ekin + \{jj)q,uj^O ■ (14) 

The expectation values are taken with respect to the BCS ground state. Note that the result 
for ps contains a contribution (Ekin, the kinetic energy per bond) which comes from a first 
order expectation value of S^. The effective action (12) obviously leads to the well-known 
Anderson-Goldstone collective mode^ for the phase excitations, with energy u{q) = v\q\, 
= Ps/k. 

One of the advantages of the present formulation is that the inclusion of the external 
electromagnetic field is straightforeward. The scalar potential leads to an extra term 

-ey(r,T)$;$^ (15) 
in the fermion action, which corresponds to the replacement 

(pj. ^ ifij. + '2ieV{r, r) (16) 
in Si^. On the other hand, the vector potential is included via the Peierls substitution 



^$:$^,exp (^--y A-d£j , (17) 
which for slowly varying A, i.e. small fields, gives rise to the replacement 

Wif ^Wip-—A , (18) 
c 

and therefore the effective action is now 

SeM = l I'' dr [ d^{K{^ + 2ieVir,r)Y + psiV^--A{r,T)Y} + Se.m. -(19) 
Z Jo J c 

Here Scm. is the action of the free electromagnetic field. One may notice that this derivation 
is entirely free of the usual "eikonal approximation" used to introduce the electromagnetic 
field. The essential point now is that (19) gives, via the Anderson-Higgs mechanism,'^ rise to 
a shift of the energy of the long-wavelength phase oscillations away from zero to the plasmon 
frequency. 

A number of comments are in order here. First, the derivation is limited to long wave- 
length, so that <^ 1. Otherwise the expectation values in (13) and (14) would have to be 
taken at nonzero q and u. Further, the derivation is limited to zero temperature, because at 
any finite temperature excited quasiparticles give rise to contributions to {nn)q^^ which are 
nonanalytic asq, uj 0. Finally, we notice that there are of course terms of higher order in (p. 
Given that V(p is proportional to the supercurrent, time-reversal invariance requires there to 
be only even powers in gradients of (p. On the other hand there is no analogous restriction for 
(p, which represents variations of the local particle density. In particular a term proportional 
to (pCVcpY is allowed by symmetry and certainly non-zero, because it represents the change 
of the kinetic energy of a supercurrent when the particle density changes. This term is in fact 
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important to obtain the London acceleration equation, as first pointed out by Abrahams and 
Tsuneto.^ 



REPULSIVE HUBBARD MODEL 

For repulsive interactions (U > 0) the Hubbard-Stratonovich decomposition in the form 
(6) is obviously unsuitable because the integrals are divergent. In this case the appropriate 
form is 



47 / dA,dAs exp \-^AAl + A^) + ieA.n + eA, 
irU J I U 



,cr. 



(20) 



Here Ac,s are real variables, and n = V'l^^ + V'J^j,' (^z = ^fV't ~ ^I^j,- the attractive 
case, one inserts this at each point in space and time and thus obtains a functional integral 
over charge and spin fields Ac,s(r, r), coupled bilinearly to the fermions. A saddle point 
approximation reproduces the Hartree-Fock results, and in particular at half-filling one 
finds an antiferromagnetic (or spin-density wave in another terminology) ground state. The 
unpleasant feature of this way of proceeding is that both Ag and A^ are scalar fields, and 
one therefore cannot construct easily the effective action for the low-energy excitations of 
the antiferromagnetic state which are spin-waves, the existence of which is of course closely 
related to the vectorial character of the order parameter. 

Alternatively, one might use a Hubbard-Stratonovich decomposition using a vector aux- 
iliary field. One then however does not even obtain the Hartree-Fock solution as a saddle 
point. A number of other, equally unsatisfactory decompositions have been discussed in the 
literature.'^ In order to obtain a spin-rotation invariant effective action for fluctuations around 
the Hartree-Fock solution, I notice*"'^ that in writing down the Hamiltonian the choice of the 
spin quantization axis is a priori arbitrary at each lattice site, and in a functional integral for- 
mulation can also vary in time. I then leave the quantization axis Qj. (r) arbitrary and integrate 
over all possible f27.(r), with the appropriate invariant and normalized integration measure 
at each point in space and time. In practice, this is achieved by introducing SU{2) rotation 
matrices Rr(r) at each point of space and time which satisfy R^. {T)azRli:{T) = Ot-(t) ■ cr. 
A convenient choice is 



R{Q.) 



e*^sin(i^) cos(i^(r,T)) 



(21) 



where 9 and 99 are the usual polar angles. I then introduce identities R^ {t)R,^{t) = 1 at 
the appropriate places in the functional integral and integrate over all configurations QriT)- 
Finally, new spinor varibles are introduced via 



$ = R+^ 



(22) 



This means that the 0-particles now have their spin along ztf2j.(T). The Hubbard interaction 
term is invariant under this transformation, and now the Hubbard-Stratonovich transformation 
can be used in its form (20) without loosing the spin excitations which are contained in the 
functional integral over Qr{T)- This also means that a nonzero saddle point value of the spin 
field As does not necessarily imply the existence of magnetic long-range order. For this to 
occur, the angular degrees of freedom have also to be ordered. 

Half-Filled Case 



At half-filling, the Hartree-Fock saddle point is antiferromagnetic for any positive U. The 
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partition function then becomes 

Z = JvnV^VS e-^«^-^"-^« . (23) 

Here Shf is the action corresponding to the saddle point, Ss represents the massive fluctua- 
tions of Ac and around their respective saddle point values, and 5^ represents the coupling 
between the angular fluctuations and the fermions and is given by 

Sn = f dT{%R^Rr<^^ - t J2 IKiRvRr' - + c.c.]} . (24) 

This term is explicitly of at least first order in time and space derivatives. An effective action 
for Qj. (r) can be obtained by integrating out the fermionic degrees of freedom, and to second 
order one finds 

'S'eff(^) = {Sn)HF — \ {S}i)HF,con ■ (25) 

Note that there are explicitly nonzero first order cumulants. Dividing Sq, into two parts Sq^^ 
and S^^r containig only time or space derivatives, respectively, one finds in particular 

{Sn,T)HF = o / moy'(-l)''(l - cos6'r)<^T- 
2 Jo ^ 

y{-\Y f AiQ) ■ dn . (26) 
^ Jo 



irriQ 



Here mo is the SDW amplitude which tends to unity for large U, Or and are the polar 
angles on site r, and A(f2) is the vector potential created by a magnetic monopole sitting at 
the center of a unit sphere.^'^ Eq.(26) can be recognized as a collection of the Berry phase 
terms for spins localized at lattice sites r. The standard procedure^'^ is now to assume at least 
local antiferromagnetic order and to write 



ilr = \ll - a^Llrir + {-lyaLr , rir ■ Lr = . (27) 

Here the antiferromagnetic order parameter field is slowly varying in space and time, and 
Lj. which represents magnetic fluctuations around g = is small, a is the lattice constant. 
Inserting this into eq.(26) one obtains (I here specialize to the two-dimensional case) 

{Sn,r)HF = ^ / drL • (n x n) . (28) 

For large U this term dominates the spin dynamics, however, for moderate or small U the 
{Sq^)hf term becomes more important. Expanding now all terms up to second order, the 
effective action is obtained as 

1 rP r 

S'eff = - / dr d^r{aiL'^ + ia2L • (n x n) + agn^ + ps(Vn)^} . (29) 

2 Jo J 

The correlation strength enters via the [/-dependence of the coefficients. In particular one 
has 

Ps = Ekin + \{j^r)q^a;^0 , (30) 

where are the transverse spin current operators. The similarity with eq.(14) for the 
superconducting case is obvious. Using the Hartree-Fock equations of motion, one can 
show that current-current correlaion functions are closely related to the spin-spin correlation 
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functions which appear in previous results. Using the notations of Schrieffer et al. and 
of Chubukov and Frenkel/^ the coefficients are given by 



Oil 
0-2 

as 



2A^ 



^-Xo+-(0,0) 



a duj 
A2 



tj=0 



xtM 



a)=0 



_ 2^^x 
a 



(31) 
(32) 

(33) 

(34) 



Here A is the SDW gap parameter, and 

1 



X 



2A3 



-E' 



sin^(A;3;a) 



4A3 



3-p2(^, |, f ; 2, 2; r) 



1 / f2 9/2 

At2^ 7?3 ^3 3-f^2l,2' 2' 2' ^' ^' 



(35) 
(36) 
(37) 



The 3F2 are generalized hypergeometric functions, and r = — 16t^/A^. 

Integrating now over L in (29) one obtains the well-known action of the nonlinear sigma 
model 



(38) 



with a = al/{Aai) + a^. From this form of the effective action we immediately obtain an 
expression for the spin-wave velocity: 



Arcryz 
A^x^ + xz 



(39) 



Numerical results for the spin-wave velocity and the spin stiffness as a function of U are 
shown in the figures. In obtaining these results it proved useful (especially for small U) to 
transform the hypergeometric functions in eqs.(35) to (37) into complete elliptic integrals, 
using formulae given by Prudnikov et al.}^ 
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FIG. 1. Spin-wave velocity (top) and spin stiffness (bottom) for the two-dimensional Hubbard 
model at half-filling as a function of U. The dashed lines are the asymptotic forms for large U, 
c = 4:'s/2t^a/U and ps = /U, respectively, which agree with lowest order results for the antiferro- 
magnetic Heisenberg model. 

One should note that the correlation strength also enters via the value of the short-distance 
(i.e. large-g) cutoff beyond which antiferromagnetic short-range order can not be defined. 
For large U this cutoff is the lattice constant, however, for general U this cutoff can be 
identified with the coherence length = t/A (the analogue of the BCS coherence length 
of superconductivity). In particular, this length becomes large for small U and diverges as 




exp{Jt/U) for U ^ 0. Consequently the coupling constant of the nonlinear sigma model 



g = clipsio) (40) 

goes to zero as [/ ^ 0, i.e. quantum spin fluctuation corrections to the Hartree-Fock 
solution become arbitrarily small for small U. The same argument also applies to the 
amplitude fluctuations contained in S^. In particular, the small-f;' limit of p., is thus expected 
to be an exact result. On the other hand, for larger U there are corrections both to c and to 
Ps, as is well-known from the large-C/ limit which is equivalent to the Heisenberg model. 
In this limit the above derivation can actually be formulated so that arbitrary (and not only 
slow long-wavelength) variations of fl^ in space and time are allowed, and one then directly 
recovers the action of the antiferromagnetic spin-1/2 Heisenberg model. ^ 

In two dimensions, thermally excited spin-wave fluctuations destroy antiferromagnetic 
long-range order at any nonzero temperature. At low temperatures the dominant fluctuations 
are always the orientational spin excitations, described by the effective action (38) (with 
possibly weakly temperature dependent coefficients). It is interesting to note that for small U 
there is a very abrupt crossover at the mean-field transition temperature between this 
regime and a high-temperature Fermi liquid state. The width of the crossover is of order 
{T^^Y /t}^ For larger U this crossover gets more and more smeared out. 

The Doped Case 

Doping away from half-filling, the antiferromagnetic ground state of the Hubbard model 
is modified quite drastically: for small U the mean-field (saddle point) approach predicts 
that for any finite doping the Neel state is replaced by a linearly polarized incommensurate 
spin-density wave.^^'^^ For small doping the magnetic structure is best described as a regular 
arrangement of linear domain walls separating commensurate regions with opposite signs 
of the order parameter. In this case the doped carriers are localized at the domain walls, 
and the systems thus remains insulating. With increasing doping, the domain wall state 
progressively transforms into a sinusoidal modulation, and this is accompanied by a metal- 
insulator transition. It would clearly be interesting to generalize the above treatment of the 
half-filled case to investigate the collective modes in the incommensurate state. 

The mean-field picture has been critized by Chubukov and Frenkel,^^ who argue that 
corrections to mean-field suppress the instability against domain wall formation at least at 
small doping, and that therefore the commensurate antiferromagnetic state survives even at 
finite doping. However, their argument is entirely based on a local stability analysis of the 
antiferromagnetic state, and this may well be insufficient when one considers domain wall 
formation which implies a global reorganization of the magnetization pattern. In fact, within 
mean-field the energy gain per doped particle due to domain wall formation is of order A,^^ 
and for small U corrections to mean-field are typically of order A^/t. This strongly suggests 
that the mean-field picture is correct for weak correlation. 

For large U an effective action for the spin degrees of freedom and the doped carriers 
can be derived, because in fact arbitrary space-time variations of f2(r, r) can be treated. For 
simplicity, one can then start from a ferromagnetic saddle point which is characterized by 
lower and upper Hubbard band separated by a gap U . For the case of electron doping, the 
chemical potential is somewhere in the upper Hubbard bands, and the lower Hubbard band 
then can be integrated out.^ In this way one obtains the effective action for the local spin 
orientation and particles in the upper Hubbard band order by order mt/JJ. To zeroth order 
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int/[/Ifind 

Seff = iyrlj2[K{dr-l^ + U)cP^-^Ml-cosA){l-cPM] 



-tJ2H^r,nr,)K<Pr'+C.C.]\ . (41) 
{rr') J 

Here refers to fermions in the upper Hubbard band, the spin index being omitted, and ^p^ , i3r 
are the polar angles of fl^. The coefficients a{Qr, ^r') come from the expression for the 
product of two R matrices: 

^^^^^ ^^^-[-e^^'a%n,-n') a%n,n') j' ^^^^ 

where 

a{nr, ^r') = lale^^""' = [(1 + ■ nr')/2]^^^ exp[zi(0^, O^,, ^)/2] , (43) 

A{Cli, r22, f^a) is the signed solid angle spanned by the vectors Jli, il2, ^3^^, and z is the 
unit vector along z. 

In the absence of particles in the upper Hubbard band, in S^^j only the purely imaginary 
term remains, which is the Berry phase of an isolated spin 1 /2, i.e., as expected, the half -filled 
Hubbard model becomes a collection of independent spins for U = 00. Introducing more 
fermions, two effects occur: (i) the factors (1 — 0*0^), previously introduced by Shankar 
from semi-phenomenological arguments^, cancel the Berry phase term whenever there is 
an extra particle on site r, i.e. one is in a spin singlet whenever two particles occupy the 
same site. Here this effect is seen directly from a microscopic calculation, (ii) the kinetic 
energy term plays a role: in particular, going around an elementary plaquette (1234) the 
lattice curl of the phases Xrr' equals $1234 = ^^2, f^a) + A(Q,:i, Q,/^, f2i)]/2, i.e. there 

is an effective magnetic field proportional to the solid angle spanned by ili, ^1,4. $1234 
is the lattice analogue of the familiar winding number density of the continuum nonlinear a 
model^^. Note that, while the gauge potential in (43) depends explicitly on z and therefore 
is not rotational invariant, the physical fluxes are. For coplanar configurations $1234 = 0, i.e. 
the phases can be removed by a gauge transformation of the ^'s. One then straightforwardly 
sees that the kinetic term is optimized by a ferromagnetic arrangement of the spins. This 
is the familiar Nagaoka phenomenon. Whether non-coplanar configurations of fi^ with a 
nonzero winding number density can lead to an energy lower than the Nagaoka state is not 
currently clear. ^° 

The first order contribution to the action is 

Slff = ^! 1% J ^ ^ (2 - <P;<P, - <PM{Qr ■ ^r' - 1) 
^ {rr') 

+ ^ [a{^lr,-^r')0i{-^r',^r")K^r" + (^■(^■] 
{rr'r") 

+ ^E(i-^:^.)^^l ■ (44) 



T" 



Here r and r" are second- or third-nearest neighbors, and the sum over r' is over all sites 
that are nearest neighbors of both r and r". In the absence of fermions in the upper Hubbard 
band, only the ^-independent part of the first term in Sl^^ contributes and represents the 
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antiferromagnetic exchange interaction between nearest neighbor sites, e.g. in this case 
Slfj + Sl^f is the action of the antiferromagnetic Heisenberg model, as already mentioned 
above. 

In the presence of a finite concentration n of extra electrons, one immediately sees the 
instability of the Neel state: in the Neel state the term proportional to t in 5°^^ does not 
contribute. If however fi^ ■ ^r' = — 1 + e^, there is an effective nearest neighbor hopping of 
order t\£\, and a corresponding gain of kinetic energy of order —t\£\n. The loss of exchange 
energy is of order e^, e.g. £ 7^ is energetically favored for any nonzero n. For an r- 
independent e one then finds a spiral, as proposed by Shraiman and Siggia^^ However, 
at least at the mean-field level discussed here the spiral state is unstable against phase 
separation.^ 

A word of caution is in place here: as describrd above, the a formalism neglects certain 
important contributions to the functional integral. This mainly comes from the fact that [ 
particle (lower Hubbard band) can equally well described by an | particle (upper Hubbard 
band), provided that the local fi^ is changed into — fi^. Ways to handle this problem will be 
discussed in a forthcoming publication. 



GENERALIZATIONS 

SU{2) X SU{2) Symmetry 

The present approach allows an interesting generalization: using the matrix 
representation^^ 

^ = i 1 (45) 

the Hamiltonian can be written in a SU{2) x SU{2) invariant way:^'^'' multiplying by 
a SU{2) matrix from the right, one generates the spin rotations discussed above {SU{2)s 
symmetry). On the other hand, multiplication from the left (SU{2)c symmetry) gener- 
ates electron-hole transformations (in the Heisenberg model this becomes a SU{2) gauge 
symmetry^^). The clear separation between the two symmetries makes the representation 
(45) very useful. 

The reason for introducing the factors (—1)'' in (45) is the identity, valid if r and r' are 
on different sublattices: 

tr[*+*,,] = + rr'si^rs) ■ (46) 

s 

Using these relations, it is straightforward to rewrite the Hamiltonian in a form invariant 
under S'[/(2)c x SU{2)s. One possible form is 

H{r, i^) = -tYl tr[*+*,,] -^J2 trk*+*,] • tr[<T*+* J . (47) 

(rr')s r 



The action becomes 

"1 

dr 

/o 



S= dr 
Jo 



(48) 



One sees immediately that 7^ breaks the SU{2)c x SU (2)s symmetry, with only U (l)c x 
SU(2)s (multiplication by a phase and spin rotation) left. 



10 



I now perform the unitary transformation = Rc'^rRs, where Rc,Rs are SU{2) 
matrices, parametrized by unit vectors fic, which vary in time and space, similarly to 
eq.(21): 

.T 



Rc (^zRc 



(49) 



A $ particle with spin up points along Clg the original (laboratory) frame. Similarly, for 
Rc ^ 1 acf) creation operator is a linear combination of the original creation and destruction 
operators. In terms of $ one has 



/ dr 
Jo 



1 



Etr[$+(a,$, + R, + <^AdrRt)Rs - ^Rc <yzR-^^r)\ 



(50) 



Now, introducing at each point in space and time an integration over flc{r, t) and il^ (r , r) 
with invariant integration measure normalized to unity, the partition function becomes 



(51) 



What has been done here is a change from a fixed reference frame in spin and particle-hole 
space to a reference frame varying in space and time. If one limits oneself to the spin rotation 
SU (2).5, this amounts to using a quantization axis varying in space and time, as before. One 
can now proceed as in the previous chapter and for example integrate out the lower Hubbard 
band. In this way one can then solve the problem mentioned at the end of the last chapter. 

The SU{2)c X SU{2)s symmetric formulation can also be used to obtain an effective 
Ginzburg-Landau type functional that contains both particle-hole and particle-particle type 
fields and recovers Hartree-Fock as its saddle point. For this purpose the Hamiltonian is 
written as 



-t Y: tr[i?,,<i?,,i?+,$,,i?+,] + {(tr[<a,$,])2 - (11^$+$,])^} (52) 

{rr') 



The non-invariant way of writing the interaction part of H is motivated by the fact that the 
saddle point equations for the resulting functional integral will reproduce the full Hartree-Fock 
equations. 

The freedom gained by the introduction of the integrations over fig and Qg will now allow 
me to restore the symmetries apparently lost in (52). First, I reintroduce the old variables ^ 
instead of Thus 



with 



dr 



(53) 



(54) 



The Hamiltonian now is 
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The SU (2)c X SU(2)s symmetry is restored due to the integrations over Q.c,s- 

At first sight, eqs. (53), (54), and (55) look like a rather complicated way to rewrite the 
initial problem. The usefulness of the new representation will become apparent introducing 
a Hubbard- Stratonovich decomposition of the four-fermion interaction. One then obtains 

Z = lvAs{r)VA,iT)Vn,{r)Vnsir)exp[-Seff{A„As)] , (56) 

where Ac,s = Ac^sQ,c,s- Note that Ac,s is integrated over (— oo, oo). The effective action for 
Ac,, is (for U >0) ' 

>5e//(A„A,) = 

r/3 



I dr^E(^c. + ^L)-ln{/^^*(^W(T)exp[-5(V'*,^,A„A,)]} ,(57) 



with 



Sir, A 



/ dr 

Jo 



^ Etr[*+(a. - /.a,)*,] - t Y: tr[*+*, 

^ 7- l^rr') 



+ W {iA,r ■ tr[*+o-*,] + A,, • tr[o-^*+*,]} 



(58) 



2 , 

The saddle point equations now reproduce the full Hartree-Fock equations: 

(dr - IJ-)lprs - t i>r's + i Yi^'^r " (T)ss'Xrs' + XK^^'' ' ^)ss''iprs' = , (59) 
r' s' s' 

|a,, = -^(tr[*>*,]) , |a,, = -^(tr[^^*+*J) , (60) 

Here Xr = (V'rs: ^{~^Ti't -s)- Obviously (59) is the effective "Schrodinger equation" for 
electrons moving in the effective space- and time-dependent fields Ac,s, whereas eqs. (60) are 
the corresponding self-consistency conditions. 

Though the Hartree-Fock equations (59) and (60) are rather standard, expanding around 
this saddle point gives rather interesting results. As far as the spin field A, is concerned, one 
finds a rather standard form. However, the transverse components of Ag contain particle- 
parrticle excitations, and one therefore obtains a functional that allows one to study the 
reciprocal effects of particle-hole and particle-particle excitations on each other. This is 
potentially useful in studying possibilities of magnetically induced superconductivity or the 
occurence or not of magnetic instabilities in relatively dilute systems, where particle-particle 
like diagrams play an important role in renormalizing the effective interactions. 

Relation with the Slave-Fermion Formalism 

The rotation matrices (21) are elements of SU{2)/U{1), and in particular have real 
diagonal elements. In the partition function (23) instead of integrating over this manifold, 
one can choose to integrate (with the properly normalized integration measure) over all SU (2) 
matrices, i.e. one writes 

Rr=( f ) , (61) 

V -Z2r 
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with the constraint |2;it. P + 1-227- P = 1- Transforming to new variables bir via 

bir = Zir{l - \4>l(}>r) 



(62) 



the constraint becomes 



(63) 



This is nothing but the constraint familiar from the slave-fermion formulation of the strongly 
correlated fermion problem: on each site there is either a spin (represented here by a boson of 
spin I {hi) or [ (62)) or a hole, represented here by a spinless fermion. Performing the same 
transformation on the action, one finds 



where the omitted terms are three-site terms analogous to those in eq.(44). This is nothing 
but the action of the t — J model in the slave-fermion representation.^'^ Note that the factors 
(1 — (/)*(/)p) and (2 — — (plicpj.,) which appeared in (41) and (44) have now disappeared. 
The same reasoning can be followed in the more general SU{2) x SU{2) formulation, and 
one then obtains other types of slave-particle boson or fermion) representations. 
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